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We derive a rigorous lower bound on the average local energy for the Ising model with 
quenched randomness. The result is that the lower bound is given by the average local energy 
calculated in the absence of all interactions other than the one under consideration. The only 
condition for this statement to hold is that the distribution function of the random interac- 
tion under consideration is symmetric. All other interactions can be arbitrarily distributed 
including non-random cases. A non-trivial fact is that any introduction of other interactions 
to the isolated case always leads to an increase of the average local energy which is opposite 
to ferromagnetic systems where the Griffiths inequality holds. Another inequality is proved 
for asymmetrically distributed interactions. The probability for the thermal average of the 
local energy to be lower than that for the isolated case takes a maximum value on the Nishi- 
mori line as a function of the temperature. In this sense the system is most stable on the 
Nishimori line. 
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O ' 1. Introduction 

Studies of random spin systems, especially spin glasses, have a long history, and a good 
amount of knowledge has been accumulated. 1,2 It has nevertheless been very difficult to es- 
tablish reliable analytical results for finite-dimensional systems, where the mean- field theory 
may not apply. A notable exception is the gauge theory of spin glasses, by which the exact 
energy can be calculated for a generic system on the Nishimori line (NL). 1 ' 3 Several other 

O ■ 

rigorous results have also been derived within the same theoretical framework. Another class 
of rigorous results relevant to the present work is a set of correlation inequalities proved for 
random spin systems under certain circumstances, which are related to the existence of ther- 
modynamic limit of various physical quantities. 4-8 It is very useful and important to find 
further rigorous relations to clarify the properties of random spin systems. 

The purpose of the present paper is to analyze the behaviour of the local energy of random 
spin systems with symmetric and asymmetric distribution functions of interactions. The first 
part of the paper concerns a generic case, not necessarily spin glasses, and therefore we do not 
use the gauge theory. We instead derive a relation to compare the average local energy with 
the same quantity calculated in the absence of all other interactions. The result is surprisingly 
simple that the former is always larger than or equal to the latter. 
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Another inequality is derived for the spin glass problem. We evaluate the probability 
that the thermal average of the local energy is lower than the same quantity in the absence 
of all other interactions. We find that this probability is a non-monotonic function of the 
temperature and assumes its maximum value on the NL. 

These two inequalities are quite non-trivial with possibilities to give us hints on the detailed 
spin states in random systems. We prove these two inequalities in the following two sections. 
The last section is devoted to discussions. 

2. Lower bound on the local energy for symmetric distribution of randomness 

Let us first define several quantities. We treat the Ising spin system described by the 
Hamiltonian 

AcV i&A 

Here V is the set of sites, and the sum over A runs over all subsets of V among which 
interactions exist. The number of sites in A is arbitrary and may be different from subset 
to subset, including the case of an external field (single-site interaction) \A\ = 1. The lattice 
structure is assumed to be reflected in the choice of A for which J a 7^ 0. 

The probability distribution of the random interaction J a is denoted by P(J A ). This 
distribution function is in general allowed to change from subset to subset, including the case 
without randomness having P(Ja) = 5{J A — J)- 

We denote the thermal average of the local energy, —J A Yli eA Si, by angular brackets as 

e A ((3,J A ) = -/jAl[S i \ , (2) 

\ ieA I p 

where (3 is the inverse temperature (3 = l/fe^T with ks being the Boltzmann constant. The 
configurational average over the distribution of randomness in interactions will be written as 

[•••]■ 

As a reference energy we define the configurational and thermal average of the local energy 
in the absence of all interactions other than the one under consideration as 

/oo 
dJ A P(J A )J A iaxih((3J A ). (3) 
-00 

The 'average' local energy will stand for the local energy averaged over both the configuration 
of randomness and thermal fluctuations. 
Our result is the following theorem. 

Theorem 1 The average local energy is larger than or equal to e A {(3) if the distribution of 
randomness is symmetric for the interaction under consideration P(J A ) = P(—J A )-' 

[e A (0,JA)]>e° A (0). (4) 
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Remark. There are no conditions whatsoever on the distribution functions of interactions 
other than Ja or the lattice structure. Thus it is possible, for example, that all other interac- 
tions are non-random and ferromagnetic. 

To prove this theorem it is useful to bring the following result into attention. 

Lemma 1 Let us divide the thermal average of the local energy into the trivial part and the 
rest: 

e A (f3, J A) = - JAtanh(/3J A ) + A(/3, J A ). (5) 

The first term on the right hand side is trivial because it is the thermal average of the local 
energy in the absence of all interactions other than J a- Then A(/3, J a) satisfies the following 
identity: 

AfR T \ J A ) , , 

J A 

Proof. Let us denote the partition function of the system for a given configuration of ran- 
domness as 

Z (P, J a) = ^2 exp I f3 Y, ^n^+^II^ I ' ( ? ) 

{S} \ BcV,B^A idB ieA ) 

where we write the interaction J a separately for later convenience. Now we evaluate the ratio 

of Z(0,-J A ) to Z(P,J A ): 

ZjP, -J a) UoflT ^ , sinh(2/3 Ja) 

= cosh(2/3J A ) + — e A (P,JA) 



1 + smh(2/3JA) ( eA (P, Ja) + Ja tanh(/?J A )) . (8) 
J A 



Similarly we obtain 



-1H {e A {/3, -Ja) + Ja tanh(/3JA)). (9) 



Z(J3,-J A ) J a 
Since eqs. ([8]) and ([9]) are reciprocal to each other, we have 

1 + sinh(^JA) 1_ 

Ja 1+ B mh(2/3J A ) 

Ja 

which leads to the desired relation ©. | 

Remark. Equation ([6]) implies that A(/3, J4) and A(/3, — Ja) have opposite signs. Note that 

the denominator of eq. ([6]) is always positive as can be verified from eq. ([8]). 

Proof of Theorem [0 We first rewrite the definition of the configurational average of 
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za(0, Ja), 

/oo 
dJ A P(J A ) [(- J A tanh(pj A ) + A(/3, J A )]' 
-oo 

/oo 
dJ A P(J A )J A ttmh(pJ A ) + [AOS, J A )} 
-oo 

= e A ((3) + [A((3,J A )], (11) 

where [•••]' stands for the configurational average over randomness of other interactions than 
J A . For the configurational average of A(/J, J A ), we have 

/oo 
dJ A P(JA)[A(/3,J A )]' 
-oo 

oo . 

dJ A (P(J A )[A((3, J A )]' + P(-J A )[A(P, -J A )]' 



/•CO 

= / dJ A P(J A )[A((3,J A ) + A((3,-J A )}', (12) 

JO 

where we have used the symmetry of the probability distribution P{J A ) = P(—J A ). According 
to Lemma 1, the sum of A(J3, J A ) and A(/3, — J A ) in the integrand of eq. (fT2j) satisfies 

sinh(2^) A(ftJ<)2 

A^ + ACft-J,)- * (13, 

1 H j Ja) 



Therefore the following relation holds, 

/•oo 

[e A (p, Ja)] = eV/3) + / dJ A P(J A ) 
Jo 



i+! 4^) A(AJj4) 

Ja 



(14) 



from which Theorem [T] follows. | 

Remark. The inequality proved in Theorem 1 is the best possible one as a generic inequality 
which is valid for any system. The reason is that a generic inequality should apply to the case 
where only the given J A exists and all other interactions are absent. We therefore have to take 
into account the characteristics of specific systems if we are to improve the lower bound. 

For example, let us consider a triangle made of three sites, between any two spins on which 
there exists the usual two-body interaction. Each of these three spins may have interactions 
with spins on other sites outside the triangle. The two-body interactions within the triangle 
are assumed to obey the same symmetric db J distribution. Then the configurational average 
of the sum of the thermal averages of these three interaction energies &tr(/3) satisfies the 
following inequality as shown in Appendix: 

[eraC/*)] > 4r(/3), (15) 
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where e^ R (/3) is the configurational average of the same energy in the absence of all other 
interactions, 

o 3 / exp(2/3J) - exp(-2/3J) , exp(2/?J) - exp(-2/3J) \ 

e T RlPJ 2 'J y exp ( 2 [3J) + 3exp(-2/3J) 3exp(2/3J) + exp(-2/?J) / ' 1 j 

This quantity satisfies the following inequality 

e%{(3) > -3Jtanh(^J) (17) 

in consistency with Theorem [TJ Equation (|15p is an improvement over Theorem [TJ which 
claims that [ctr(/3)] > — 3Jtanh(/3J), for the special case of the triangular arrangement of 
three sites. 

3. Local energy for asymmetric distribution of randomness 

We next analyze the probability distribution of the local energy for asymmetric distribu- 
tion function of random interactions. 

Assume that all the interactions in the system share the same distribution function satis- 
fying 

P(-Ja) = exp(-2/3 N J A )P(J A ), (18) 

where /?n is a parameter corresponding to the inverse temperature on the NL as will be 
exemplified below. This constraint is known to allow us to apply the gauge theory. 1 ' 3 For 
example, in the case of the ± J model 

P(J A ) = P 5(J A - J) + (1 - p)5(J A + J) (j>0,p>^\, (19) 

we find 

1 — p 



exp(2/? N J) = (20) 



For the Gaussian model 



P(JA) = ^=ex P {- iJ \ J 0) \ (./„>(.,. (21) 



Pn is 

R,r — 

J 2 



/3n = 4 (22) 



Our result is the following theorem. 

Theorem 2 The system defined above satisfies the following inequality, 

r hp, j a) 

[\ A(/3,J A ) 
with A(/3, Ja) defined in eq. |5]). 



< 



A(/3 N , J A ) 



(23) 



Remark. The ratio 

AQfl, Ja) 
I A(/3, J A ) 



(24) 
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is +1 if the thermal average of the local energy ca(P, J a) is lower than the trivial value 
— J a tanh(/3J 7 4) and is —1 otherwise. Thus the configurational average of this ratio is the 
difference of p\ and P2, where p\ is the probability that ca(/3, J a) is lower than — J a tanh(f3 Ja) 
and P2 is for the other case. Since P1+P2 = 1, the difference satisfies p\ — p2 = 2p\ — 1. Hence 
the configurational average of the above ratio is a measure of the probability p\ that e J 4(/3, J a) 
is lower than — Ja tanh((3 J a)- 
Proof. We rewrite the ratio as follows: 



AQ3, J A ) 
A(/3, J A ) 



oo 
oo 



/oo 
P(Ja 
-oo 

p(-Ja) 



A(A Ja) 
A(A J A ) 

" A(/?,-Ja) 
.1 A(^-Ja) 

P(J A )exp(-2/3 N J A ) 

A(/?,-Ja) 



AQS, - J A ) 
A(A-Ja) 



Equation ([6]) implies 



exp(-2/3ArJA) 
A(/3,-Ja) 



I A(/3,-J A ) 
A(/3, Ja) 



Thus, we obtain 



A(/3,Ja) 
ACS, Ja) 



a(A-Ja)| 

exp(-2/3 N JA) 



A(/3, Ja) I 
-A(/3,Ja) 



(25) 



(26) 



AGS, Ja) 
cosh(2/3 N Ja) - sinh(2/? N Ja)) 



-A(/3,Ja) 
A(/3, Ja) I 



According to the gauge theory, 1,3 we may perform the local gauge transformation 

Si ► SiO~i Ja 



(27) 
(28) 



J A \\ 0~i, 

where <jj is a local gauge variable which takes the values ±1 (arbitrarily fixed at each site). It 
is useful to note that A(/3, Ja) is invariant under the above gauge transformation. Then, we 
obtain 



cosh(2/3 N J a) - sinh(2/3 N J A ) \\ a% 



A(/3, J a) 
A(AJa) 

cosh(2/3 N J A ) - sinh(2/? N Ja) ( \\ o 
sinh(2/?N Ja) 



-A(AJa) 



/3n/ 



-A(/?,Ja) 
I AOS, J A ) 



A(A Ja) 



1 + 



Ja 



-Ja m 



O + Ja tanh(/3 N Ja) 

i^A I p N 



-A(/3, Ja) 
A(A Ja) I 



Hence we have the following identity, 



r a(aja) " 




L|A(/3,Ja)|. 





^A 



A(A Ja) 



A(/3, Ja) 



(29) 



(30) 
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from which we find 

AOS, J a) 



A(A J a) 
We therefore have 



A(/3,J A ) 



A(/3,Ja) I J 



sinh(2/3 N J A ) A(/3, J A ) 

"A(/} N , J^j: 



-/.4 



A(/3, Ja) 



A(AJa) 
A(/3, Ja) 



sinh(2/? N Ja) am r s A(/?,Ja) 



2Ja 



A(p\ Ja) 



When f3 = (corresponding to the NL), we get 



A(/3n,Ja) 
A(/? N , Ja) 



sinh(2/3N Ja) 
2J^ 



A(/3 N) Ja) 



(31) 



(32) 



(33) 



On the other hand, for an arbitrary temperature, we obtain the following inequality from eq. 



AQg, Ja) 
AOS, Ja) 



< 



sinh(2/3 N Ja) A(/3,J a ) 
2Ja A(/3n ' Ja) |A(/?,Ja)| 

sinh(2/?N Ja) 



2Ja 



A(/3 n , Ja) 



(34) 



From eqs. (|33|) and (|34|) . we derive eq. 
Remark. The quantity — Ja tanh(/3 Ja) is the thermal average of the local energy in the ab- 
sence of all other interactions (isolated case). Theorem [2] states that the deviation A(/3, Ja) 
from the isolated case is more likely to be negative on the NL than at any other temperatures. 
In this sense the system is most stable on the NL than at any other temperatures. This obser- 
vation reminds us of a similar non-monotonic behaviour of the spin orientation as discussed 
previously. 9 

4. Summary and discussion 

We have proved two inequalities on the behaviour of the local energy for the Ising model 
with random interactions. The first inequality concerns the case of symmetric distribution 
function of randomness and states that the configurational and thermal average of the local 
energy is not smaller than that of the same local energy in the absence of all other interactions. 
Therefore an introduction of any interactions to the isolated case always increases the average 
local energy. This is to be contrasted to the result for ferromagnetic systems, the Griffiths 
inequality, 10 in which an introduction of other interactions always lowers the thermal average 
of the local energy. The random system under consideration behaves quite differently from the 
conventional ferromagnetic systems. Notice that our result is a special case of the conjecture 
presented by Contucci and Lebowitz. 8 

The second inequality is for asymmetric distribution of random interactions. The prob- 
ability of decrease of the thermal average of the local energy from the isolated case is a 
non-monotonic function of the temperature and assumes the largest value on the NL. This is 
again a somewhat counter-intuitive result. Note that the average of A(/3, Ja) itself, not the 
sign of it, vanishes on the NL. 1,3 The sign of A(/3, Ja) may be positive or negative depending 
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on the configuration of randomness. If we focus ourselves only on the sign of A(/3, Ja), ig- 
noring the magnitude, we find that the probability for this quantity to be negative is largest 
on the NL. The NL therefore occupies a special position in the phase diagram from the view 
point of spin configurations. Further clarifications are necessary. 
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Appendix: Derivation of eq. (j!5|) 

In this appendix we explain the derivation of eq. (|15p . Let us consider the energy of a 
three-spin system consisting of spins S\, S2 and S3, which interact with each other by three 
two-body interactions. We number the configurations of three interactions, {Ji2> ^23) ^31}) 
namely, {J, J, J}, {J, - J, - J}, {-J, J, - J}, {-J, - J, J}, {-J, -J, -J}, {-J, J, J}, {J, - J, J} 
and {J, J, —J} as configuration 1 to configuration 8. Note that configurations 1 to 4 are 
unfrustrated while 5 to 8 are frustrated. We denote the partition function for configuration n 
as Z(/3,n) and the thermal average of the energy of a bond {ij} as ejj(/?,n). Note that the 
interactions with spins outside the triangle are arbitrarily fixed. 

It is convenient to evaluate the ratio of two configurations: 



Z((3,l) 



((cosh(2/3J) - sinh(2/?J)S , 2 S 3 )(cosh(2/3J) - sinh(2/3J)S 3 Si)) pi 
cosh 2 (2/3J) + cosh(2/3J) sinh(2/3J)^^ ' ' ^ ' ' 



J 

_ sinh 2 (2/ 3J)£^M, (Al) 

where (• ■ - )bi denotes the thermal average at inverse temperature (5 under the bond config- 
uration 1. 

Now we note the trivial inequality, 

Substituting the results obtained by repeating the procedure to derive eq. (|A-ip into eq. (jA2]), 
we find 

t=l i=l 

where 

eTR((3,i) = e 12 (/3,i) + e 23 (A0 +e 3 i(P,i). (A-4) 
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From eq. (|A-3j) . it follows that 

if (R -w 7 exp(2/?J) - exp(-2/3J) 

12 ^ 6TR(/? ' * } " ~ J exp(2/?J) + 3exp(-2/?J) ' (A " 5) 

where the right-hand side is equal to the energy of a non-frustrated triangular cluster in 
isolation. Similarly, it holds that 

Z(/?,6) 



■^-^ = ((cosh(2/3J) + sinh(2/?J)5 2 S 3 )(cosh(2/3J) + sinh(2/3J) < S 3 5i)) /3j5 
= cosh 2 (2/3 J) + cosh(2/3J) sinh(2/3J)- 



2/o«* , „^o«T^ ; ^oAn e 23(/?,5)+e 3 i(/3,5) 



J 

+ sinh 2 (2/3J) ei2( f 5) , (A-6) 

where (• • ■ )/? 5 denotes the thermal average at inverse temperature j3 in the bond configuration 
5. By the same procedure as above, we obtain 

8 8 

12 cosh 2 (2/3J) + 2 cosh(2/3J) sinh(2/3 J) ^ fl^ill + sin h 2 (2/3J) ^ fl^ill > 12 . (A-7) 
Consequently, 

8 



If m .s^ T exp(2/3J) - exp(-2/?J) 



12 v ' ' ~ 3exp(2/3J) + exp(-2/3J) 

where the right-hand side is equal to the energy of a frustrated triangular cluster in isolation. 
Finally, we obtain that 

i=l 

3 / exp(2/?J) - exp(-2/3J) exp(2/3J) - exp(-2/3J) \ 

2 V ex P(W) +3exp(-2/3J) 3exp(2/3J) +exp(-2/3J)y ' 1 ' } 



This is eq. (115]) . 
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